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Abstract 

The majority of literature on robust control assumes 
that a design model is available and that the uncertainty 
model bounds the actual variations about the nominal 
model. However, methods for generating accurate 
design models have not received as much attention in 
the literature. The influence of the level of accuracy of 
the uncertainty model on closed loop performance has 
received even less attention. The research reported 
herein is an initial step in applying and extending the 
concept of model validation to the problem of 
obtaining practical uncertainty models for robust 
control analysis and design applications. An extension 
of model validation called 'sequential validation' is 
presented and applied to a simple spring-mass-damper 
system to establish the feasibility of the approach and 
demonstrate the benefits of the new developments. 

Introduction 

Robust control theory guarantees that a feedback 
control system can be designed that will maintain 
desired levels of stability and performance subject to 
modeling errors and uncertainties. 11 '" 1 There is, 
however, an underlying assumption that the 
uncertainty model used in the design effectively 
characterizes the differences between the responses of 
the true system and the nominal design model. This 
means that the family of responses associated with the 
design model contain the responses of the true system. 

It is impossible to conclusively prove that this 
assumption is satisfied for any real system. However, 
if there is sufficient knowledge of the response 
characteristics of the true system then a model can be 
generated that, subject to the available knowledge, 
characterizes the possible range of responses that can 
be produced by the system. The goal is to 
systematically generate such a model with the added 
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property that it be accomplished with minimal 
conservatism. That is, the model should characterize 
responses of the real system that are possible but 
should not characterize responses that are not 
possible. Of course, it is impossible to conclusively 
prove this as well. However, if the available 
knowledge is comprehensive enough it is possible to 
achieve this goal within some qualitative level of 
confidence. 

Model Validation 

A concept called "model validation" has recently 
been developed to attack this problem. The idea 
behind model validation is that given input and output 
data for a system that is otherwise unknown, a model 
can be generated that driven by the same input (and 
possibly some additional inputs) can exactly reproduce 
the output of the true system. The frequency domain 
version of this statement is depicted in block diagram 
form in Figure 1 . 

The input to the true system u(s) produces the 
output of the true system y(s) . The same input drives 
the model along with an external disturbance w(s) . 
Another external input v(s) is added to the response of 
the model to produce the validation output y (s) . The 
external inputs are included because the true system 
output includes artifacts of the method(s) by which it 
is obtained including external disturbances, estimation 
errors, and quantization effects. The external inputs 
provide a mechanism to account for these effects. 
Validation is achieved by choosing the model and the 
external inputs w(s) and v(s) so that the difference 
between the validation output and the true system 
output e(s) is identically zero for the available 
input/output data. (The Laplace variable S will be 
omitted henceforth for ease of discussion. ) 



Figure 1 - Model validation block diagram. 
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The goal is to choose the validating model and 
external inputs in such a way that the resulting model 
can be used for robust control design. There are 
infinitely many ways to choose a validating model and 
so the next step is to develop a means for generating a 
validating model that can be used in robust control 
design and demonstrate reduced levels of conservatism. 
Robust Control Models 


Robust control design requires a linear design model 
that characterizes parameter variations, unmodeled 
dynamics, and/or other uncertainties relative to a 
nominal model. The common format for describing 
such a model in block diagram form is shown in 
Figure 2. |L ~ The nominal model is represented by the 
(2,2) element of the transfer matrix P while the other 
elements represent various aspects of the uncertainty. 
The matrices £2 and A g characterize the structure of 
the uncertainty in the model; £2 is a real diagonal 
positive definite weighting matrix and Ag is a block 

diagonal matrix with unity norm (a so-called unit ball). 
The elements of A g are only specified as to whether 

they are real or complex valued, satisfy a particular 
block diagonal structure, and satisfy the norm bound 
condition. 

For a given set of inputs U, V, and W, the output of 
the uncertain model y represents a family of outputs 
determined by all possible choices for Ag that satisfy 
the prescribed structure. One measure of conservatism 
for the uncertain model described above is the size of 
the external inputs V, W and the norm of uncertainty 
weighting matrix £2. A reasonable goal would be to 
generate a validating model that minimizes these 
quantities. 

The model validation problem can be broken down 

into several sub-problems. First input/output data from 

the true system must be obtained. It may come from 

experiment, actual operation, simulation, or computed 

directly from a high-fidelity model that is used to 

represent the true system. The second sub-problem is 

choosing a structure for the model. Robust control 

design relies on a linear model called a linear fractional 
n 21 

transformation (LFT) ' and is represented by the 
block diagram in Figure 2. The equation for the LFT 
in Figure 2 is 


P 23 + P 21 A B Q [ l_P 11 A B a ] 1p 13J^ 


p 22 + p 2lAg£2[l- p iiAgQ] Pi 2 


w 

IU 


( 1 ) 



u(s) 

v(s), w(s) 


Figure 2 - Robust control model block diagram. 

Choosing the model structure involves selecting the 
F^j , the block structure of £2 and A g , and the external 

inputs V and W. The P|j arise from the form of the 

nominal model and the representation of uncertainty 
(e.g., real or complex valued parameter uncertainty, 
additive uncertainty, input or output multiplicative 
uncertainty in unstructured or structured forms). 
Subsequently the form of £2 and Ag are chosen 

consistent with the uncertainty representations. The 
third sub-problem is setting up an optimization 
problem that solves for the smallest external inputs 
and uncertainty weights that, subject to the choice of 
model structure, makes the validation error e (from 
Figure 1) zero. 

Extending the Model Validation Method 

The current formulation of the model validation 
problem optimizes on a single scalar uncertainty 
weight with a bound constraint on the magnitude of 
the external inputs. 131 There may be other formulations 
that may be better able to incorporate the available 
knowledge. An alternate, more general, approach is 
proposed and developed in this paper. 

The choice for the structure of the uncertainty plays 
a major role in the relative conservatism of the 
resulting validating model. How does one choose the 
uncertainty structure upon which the method operates 
to get the "best" model? It may be possible to 
parameterize the potential uncertainty model structures 
and utilize optimization to select the best structure. An 
approach to address this issue by means of sequential 
uncertainty modeling is proposed and developed here as 
well. 

A brief review of the key technical aspects of the 
model validation method will be presented first. This 
will be followed by a discussion of the extensions that 
have been made to the method as part of the current 
research. This is followed by simple examples that 
demonstrate the application of the model validation 
method to generate an uncertainty model for a simple 
dynamic system. 
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The perspective on model validation and uncertainty 
modeling addressed here is a direct extension of the 
work of Lim and Giesy. 131 The model validation 
problem statement and solution is theirs and the 
current work builds on that foundation. It is 
summarized and interpreted here to provide a basis for 
the extensions described subsequently. 

Model Validation Problem 13 ' 41 

Assume that input/output data for a physical 
dynamic system is available and that a robust control 
design model for the system is desired. Assume further 
that a linear model, specified uncertainty structure, and 
a specified level of external disturbance is available for 
the same system. The model elements can be described 
as a linear fractional transformation (LFT) as shown in 
Figure 2. The input data for the physical system can be 
used as input to the LFT to produce another set of 
output data. Model validation is achieved by selecting 
certain elements of the LFT and its external inputs in 
such as way as to make the difference between the 
outputs of the physical system and LFT model exactly 
zero at each frequency. Figure 1 depicts the model 
validation problem in block diagram form. In equation 
form model validation is simply 

0 = e= y-y =e° -P 23 j^} -P 2 l^ ( 2 ) 

where 

ey = y - P 22 U 


Recall that the nominal linear model P 2 2 and the 
input data U are given. Therefore the quantities 

It t1 t * 

P231V W r and P21^ can be used to 


characterize the differences (errors, uncertainties, etc.) 
between the true system outputs and the model 
outputs. Validation is achieved by solving Eqn. (2) for 
these quantities subject to the specified structure of the 
model. 

The model structure is described by the LFT in 
Eqn. (1). The specified uncertainty structure is 
represented by F^j , i,j =£ 2 and A g . The remaining free 

variables are V and W and the uncertainty weights in 
G . Note that G has a diagonal structure dictated by the 
structure of Ag which is block diagonal. There is one 
free weight C 0 j associated with each block of Ag as 
shown in Eqn. (3) where the l n ,n = 1 ,2,..., x are 
identity matrices with dimension equal to the 


dimension of the Agj , i = 1,2,..., x and X is the 
number of uncertainty blocks. 


G = 


co-il 


1'1 


A B = 


CO21 2 


®x^x J 


A B 1 


Ab, 


Ab 


(3) 


Each Agj are unit norm matrices with specified 

structure (diagonal or full, real or complex) but whose 
elements are otherwise arbitrary. Therefore the free 
variables that can be selected to achieve model 
validation are V, W, ft)j,and Agj , i = 1,2,..., x. 

The current model validation problem formulation 
combines V and W into a single vector. In addition, it 
is assumed that these signals satisfy the Euclidean 
norm ' bound constraint 



This comes from the notion that the external inputs 
are associated with some random phenomena that are 
unknown except that they are stable and bounded. 

Model Validation Solution 131 

The solution of the model validation problem 
described above is accomplished using a linear vector 
space approach. The model validation equation, 
Eqn. (2), can be rewritten 



where M = [P 21 P 23 ] and ey , P 21 . and P23 are 
known. A necessary condition for Eqn. (5) to be 
satisfied is that ey lie in the subspace spanned by the 
columns of M , that is 

ey elm(M) (6) 

The general solution of Eqn. (5) can be written as 


Note that the superscript T is used herein to denote 
complex conjugate transpose. 


' All norms referred to herein are Euclidean norms. 
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% 

V 


[= M + ey +N M e 


wj 


where 



coldim(M) >rowdim(M) 
col dim(M)< rowdim(M) 


(7) 


and where the columns of Nm form a basis for the 

null space of M . Note that the vector 0 is arbitrary 
because when it is premultiplied by Nm the product 

lies in the null space of M which when premultiplied 
by M , as in Eqn. (5), is identically zero. However, 0 


does contribute to the norm of 


w 



. Also 


recall that the vector formed by V and W has a norm 
bound constraint. Therefore, 0 is arbitrary only to 
the degree to which the norm bound is satisfied. 

In order to assure that this constraint is satisfied 
consider a partition of Eqn. (7) 



+ (Nm ) e 

VW 


( 8 ) 


where the subscript VW represents taking the rows of 
the respective matrix associated with V and W. Now 
(Nm) can be represented in terms of its singular 


value decomposition by 


(n m ) =[Ui u 2 


fSl 0] 

1 

_< 
— 1 

o' 

O 

_v 2 t J 


(9) 


The columns of the matrices V-| and V 2 span the 

image and null spaces of (Nm) , respectively. 

vw 

Therefore, 0 can be written as the sum of two 

vectors, one in the image of (Nm ) and one in the 

vw 

null space of (Nm) . 

vw 


0 = Viy + V2\|/ 


( 10 ) 


Eqn. (8). The properties of the singular value 

decomposition of (Nm) allow considerable 

vw 

simplification and the elimination of \| / from the 
expression. 

W = (m + ) eJ+UiSiY (11) 

(WJ ' 'vw 


Recall the norm bound on V and W from Eqn. (4) 
and note that the matrix [Ui U 2 ] is unitary so that 

the constraint in Eqn. (4) can be written 


i 

CM 

> 

1 

1 

j v j 

|w| 

KJ 

[wj 


( 12 ) 


Premultiplying Eqn. (11) by [Ui U 2 f, 
substituting it into Eqn. (12) and simplifying the 
expression results in 


. 

2 


K(M + ) e°+SiY 
| v V w 

<1- 

M M+ L e ? 


(13) 


This condition can only be satisfied when the right 
hand side is non-negative which further implies that 



(14) 


This condition and that imposed by Eqn. (6) are 
required for the existence of a validating model. If these 
conditions are not satisfied then the given nominal 
model and uncertainty structure cannot validate the 
input/output data. In such cases the nominal model 
and/or uncertainty structure and/or the definition of the 
external inputs must be modified. 

Now define a variable (f> such that 

4 = u|(m + ) eJ + Siy. (15) 

The condition in Eqn. (13) then becomes a constraint 
on (f> such that 


Note that because \|/ multiplies V 2 which spans 
the null space of (Nm) it does not contribute to the 


norm of 


w 


and is therefore completely arbitrary. 


However, y does contribute to the norm of 


and 


is therefore constrained. To determine the constraint on 
y substitute the expression for 0 from Eqn. (10) and 

the expression for (Nm) from Eqn. (9) into 

vw 


M 2 


<1 



Solving Eqn. (15) for y gives 

r = S^X(M + ) vw e°]. 


(16) 


(17) 


There is an additional constraint on the selection of 
(f> and \|/ associated with the structure of the 
uncertainty described in Eqn. (3). Assume that tf) and 
\|/ satisfy the necessary conditions for existence of a 
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validating model. Eqn. (7) can then be written in terms 
of (|) and \|/ as follows. 

= [wi + -N M V 1 S 1 -1U 1 T(M + ) vw 1 je5 

(18) 

n m v 2 ]|*} 


q 

• v 
w 


Note, however, that the vector q is not an external 
input, it depends directly on the uncertainty structure 
(see Figure 2) and needs to satisfy the following 
relations. 


n=[Pn 


c, = Aeflr| 



(19) 

( 20 ) 


Substituting Eqn. (18) into Eqn. (20) results in the 
following equation. 


n= p*n 

+ Pl 2 u + [Pii P^lfNM^Sr 1 N m V 2 


( 21 ) 


Therefore, given (() and \|/, the model validating 
solution is computed by Eqns. (18) and (21) subject 
to the further constraint that Eqn. (19) be satisfied. 

Satisfying Eqn. (19) assures that the validating 
inputs determined by (|) and \|/ , U , and ey also 
generate Tj and c, that satisfy the given uncertainty 
structure. Once such a solution is determined an actual 
uncertainty A g that achieves validation for the given 
input/output data and system model can be computed 
using Eqn. (19). Note that this computation is not 
required during the optimization used to compute ([) 
and \|/ . In addition, the solution to this equation is 

usually underspecified. That is, there are generally 
many possible choices for the elements of Ag that 

satisfy the equation. A method for computing an 
optimal validating uncertainty will be presented 
subsequently. 

In summary, the Lim-Giesy solution of the model 
validation problem boils down to (1) obtaining 
input/output data for the true system, (2) selecting a 
nominal linear model, a particular uncertainty 
structure, and set of external inputs, (3) testing whether 
the candidate model structure can generate a validating 
set of inputs and uncertainties, and (4) choosing (|) and 
\|/ such that the model is validating, the external 


inputs are appropriately bounded, and the uncertainty 
structure is not violated. The following sections will 
build on this solution approach. 


Optimization Problem 

The largely arbitrary nature of (|) and \|/ suggest 
that constrained optimization can be used to select 
them subject to the necessary constraints and additional 
user specified constraints and objectives. The necessary 
constraints are the bound on ([) and the requirement 

that the uncertainty structure be satisfied. Additional 
constraints and objectives may be chosen to minimize 
the magnitude of some or all of the external inputs, 
minimize the magnitude of the uncertainties required to 
achieve validation, or trade off the magnitude of 
particular elements of the uncertainty structure. 

The approach described in the literature 13 ^ is to place 
a scalar weighting factor p on the uncertainty 
weighting matrix G= pGo (with Go specified) that 
is minimized by selection of 4> and \|/. In this way, 
the magnitude of the contribution of the uncertainty 
model to the validating response is minimized. This 
approach almost guarantees that the external inputs 
will take on values such that the norm bound condition 
in Eqn. (4) is active. This approach is fine if the 
magnitude of the external disturbances are known or if 
the uncertainty model is very simple (e.g., a single 
block). However, this is not usually the case. 

A more comprehensive approach is to optimize over 
the weights on the uncertainty blocks coj with the 
added objective of minimizing the magnitude of some 
or all of the elements of the external inputs V and W. 
The cost function may be defined so as to seek the 
smallest possible norm of the uncertainty weight 
matrix G and of tf) (which is essentially equivalent to 

minimizing the norm of the external inputs V and w) . 
This is the approach that has been implemented in the 
current research. 

In more precise terms, assuming that none of the 
uncertainty blocks are associated with repeated scalar 
uncertainties, the problem statement is as follows. At 
each frequency for which input-output data is available 

min 

4b v. ®j 


®i 


©1 


IIq 2 


( 22 ) 


subject to the constraints 



(23a) 

(23b) 
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where Qi and C>2 are diagonal semi-definite 
weighting matrices and 10 Iq is the norm of the 
vector (•) with respect to the weighting matrix Q . 

At each frequency the objective function seeks to 
minimize, in a normed sense, the contributions of the 
various blocks of the uncertainty model while 
simultaneously minimizing, again in a normed sense, 
the size of the external inputs subject to two sets of 
constraints. The first constraint Eqn. (23a) insures that 
the validating inputs t, , V and W are consistent with 
the structure of the uncertainty model determined by 
the ©j andABj in Eqn. (3). That is, there exists an 

admissible uncertainty satisfying the norm bounds and 
block diagonal structure imposed by the uncertainty 
model. The second constraint Eqn. (23b) insures that 
the validating external inputs satisfy the bound 
condition from Eqn. (4). 

If the uncertainty model further specifies that the 
structure of some of the uncertainty blocks represent 
repeated scalars such that 

ABj = Sfn; (24) 

where l nj is the identity matrix of order rij, then 

additional constraints must be imposed. The 

additional constraints are equality constraints and can 
be written 

^j=©j8itii (25) 

for those blocks that represent repeated scalar 

uncertainties where qj and T)j are vectors 
corresponding to the ith uncertainty block A Bj • The 
scalar 8j can be real or complex valued depending on 
the uncertainty description but are always norm 
bounded by unity consistent with the definition of 

a b . 

Sequential Model Validation 

A further extension of the model validation approach 
is possible if elements of the uncertainty model is 
specified (i.e., one or more blocks of A B ). The effect 
of the specified elements on the model validating 
solution can be identified and eliminated from the 
solution of the remaining validating signals. This is 
done using an image space-null space decomposition 
similar to that used in the derivation of the model 
validation solution. 

Consider the situation in which the uncertainty 
structure is partitioned into two parts, one of which is 
assumed to be known a priori, f2|<AB k > an d one that 

is to be determined, £2 U A Bu • The constraint in 
Eqn. (19) can be written as two constraints on the 


elements of £, and q associated with the two 
uncertainty model components, the first known and the 
second unknown. 

= As k Hk =A k % (26) 

= ABu Hu = Au Hu (27) 


The fact that f^kAB k is specified influences the 
possible choices for ([) and \|/ to achieve validation. 

Consider the partitions of Eqns. (18) and (21) 
corresponding to the known uncertainty block. These 
can be written in terms of a linear equation in ([) and 

¥ ■ 


^ k = ^+R k 


'lk - hk +^k 


where 




(M + ) k -(N M ) k V 1 Sr 1 Ui T (l 

Hk = ( p 12)k u + [( p 1l)k ( P 1 3 ) k 

m + -n m v 1 s 1 - 1 u 1 t (m+) 


(28) 


M 


l e y 

vwJ 


and where 


Rk = (N M ) k [VlSr 1 V 2 ] 

T k = [(Pll) k (Pl3) k ] N M[ViSr 1 V 2 ]' 

The notation (•) k represents the rows of (•) 

corresponding to the known uncertainty block and its 
inputs % and outputs q k ■ 

Substituting Eqn. (28) into the constraint in 
Eqn. (26) and solving for ([) and \|/ results in the 
following expression. 

{*} = x k(^k~ 4 l<’lk) ™ 

where 

X k = AkTk -Rk 


X k ' 


Xk T (x k X k T 

(x k T x, 


-1 


x k T . 


coldim(Xk) > row dim(Xk ) 
coldim(X k )<row dim(X k ) 


The right hand side of Eqn. (29) can be also written 
in terms of two components, one that lies in the image 
of X k and one that lies in the null space of X k . The 
singular value decomposition of X k provides a means 
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by which bases for these components can be 
determined. 


Xk=U Xl Ux 2 


( S X1 °1 

L 0 oj L v X2 Tj 


V Xl T l 


(30) 


The columns of the matrices v Xi and Vx 2 span the 
image and null spaces of X k , respectively. Therefore, 
an expression for (|) and \|/ can be written 



V Xl a+ V X2 P 


(31) 


It is readily shown that substituting Eqn. (31) into 
Eqn. (29) simplifies to the following. 

a = S Xl _1 U Xl T (^-A k q°) (32) 

Note that the solution is independent of P . Also note 
that substituting Eqn. (32) into Eqn. (31) results in 
the following expression for ([) and \|/ in terms of P . 

{J}= X ;(^- i k'lS) +V X2P (33) 


As a result, the constraint in Eqn. (26) is satisfied 
by Eqn. (33) and P can be chosen to satisfy the 
remaining model validation relations and constraints. 
This is accomplished by substituting Eqn. (33) into 
the partitions of Eqns. (18) and (21) associated with the 
uncertainty block! s) remaining to be determined. The 
new model validation problem becomes finding a 
solution for P that satisfies the combination of the 
following equations and the constraint in Eqn. (27). 


c u 


[KL-< N »>U v ' s f 1u,T K 


+(NmUw[W V 2 ]x k + (^-A kllk °) LU ’ 

+ ( N M \j V w [ V 1 S 1 _1 V 2Fx 2 P 


same computing machinery used to solve the original 
model validation formulation can be used. 

This new result provides a means for several 
variations of model validation that can be used in the 
solution of practical issues in uncertainty modeling. 
For example, the selection of the appropriate level of 
parameter uncertainty can be chosen by comparing (or 
optimizing) the level of additional uncertainty needed 
to achieve validation. This concept is developed further 
in the example section. 

Minimum Norm Validating Uncertainty 

The model validation based uncertainty modeling 
method provides uncertainty weights. It is useful to be 
able to compute an actual validating uncertainty (i.e., 

and Ab ) that has the same weights. 

There are basically two types of validating 
uncertainty models to consider - one associated with 
repeated scalar blocks and one associated with full 
complex blocks. Both types result from the solution of 
the constraints in Eqn. (19) for a specific £2 Ab given 

an input-output pair and validating uncertainty 
weights. Recall that the validating uncertainty model is 
determined by solving for ([) and \)/ . The values of q 
and c, can then be computed using Eqns. (18) and (20) 
and used to solve Eqn. (19). 

Solving for a validating uncertainty associated with 
a repeated scalar block is trivial because the number of 
unknowns is the same as the number of equations, as 
is evident in Eqn. (25). 

The solution of a validating uncertainty for a full 
complex block is more complicated because there are 
typically more unknowns than equations. One way to 
obtain a solution for a full complex block is via 
optimization. A family of validating uncertainties can 
be parameterized and these parameters can be optimized 
to satisfy Eqn. (19) and match the desired uncertainty 
weight. 

Consider a set of candidate validating uncertainty 
blocks 


1 lu=( p i2) uvw u 

+ [( p iiWw (P^UwlK-NM^SC^^irjJeO (35) 
+[( p iiUw (PiaUw^ViSr 1 V 2]Xk + (^k-Akh k °) 
+ [( p n) U vw (^UwM^sr 1 V2 Kp 


a'= 0 

. §' ... o 

, i = 1,2,...,n 

(36) 

where 

5li] 

s' =■ ; 

Sni 1 

1 

(37) 


Notice that the form of the expressions in Eqns. 
(34) and (35) is very similar to those in Eqns. (18) and 
(21). Each expression is the sum of a constant term 
and a term multiplying an unknown vector ( c|> and \|/ 
in the earlier case and [3 in this case). As a result, the 


Substituting each block into Eqn. (19) simplifies to 
the following expression for the 8 1 . 

8'=^., i= 1,2, ,n (38) 
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The sum of all the candidate solutions A 1 scaled by 
constant factors kj results in the following equation. 

^nr-L r ki{ §1 } - k n{ §n }]n (39) 

£ki 
i= 1 

Equation (39) describes a family of validating 
uncertainty blocks for arbitrary choices of the kj . 

A minimum norm validating uncertainty can be 
obtained by choosing the kj to minimize the difference 
between the specified uncertainty weight CO and the 
norm of the uncertainty matrix described in Eqn. (39), 
or equivalently 

min. ,~,i n~ 

- | co— | A| | subject to £kj=1 (40) 

K i i=1 

where 

A=£kjA' (41) 

i=1 



i 


b=J 


C 2 >k; 


m 2 



1=1 


C 3 >k; 


Figure 3 - Mass-spring-damper system. 


A is a minimum norm validating uncertainty. 

Examples 

Two examples will be presented to (1) demonstrate 
the result of optimizing the weighting for each block 
within a specified uncertainty structure, and (2) 
demonstrate the application of sequential model 
validation when part of the uncertainty model is known 
or otherwise specified. 

System Model 1 351 

Consider the simple physical system depicted in 
Figure 3. Two masses are connected to ground and each 
other via springs and dampers. The values of the 
masses, spring stiffness, and damping ratios are 
presented in Table 1. The inputs to this system are 
force actuators on the two masses and the outputs are 
the velocities of the two masses. 

The equations of motion for the system in Figure 3 
can be written 


mi 

0 


0 

m2 


xi 
X2 
"k-j + k 2 
_ - k 2 


X 1 

x 2 


Oj + C2 — C2 

c 2 C2+C3J 


k 2 + k 3 J[x 2 J [f 2 


xi 

X2 

h 


(42) 


These equations of motion can be written in modal 
coordinates and transformed into the frequency domain 
as a basis for demonstrating the model validation 
solutions described in the previous sections. The full 


Table 1 - Parameter values for mass-spring-damper system. 


Mass (kg) 

Damping Ratio 

Stiffness (N/m) 

m. 

m 2 

Ci 

C 2 

C 3 

k. 

k 2 

^3 

0.45 

0.45 

0.40 

0.40 

0.40 

150 

150 

150 


order system will be used as the "true system" shown 
in Figure 1 . The nominal linear model is obtained by 
truncating the second (higher frequency) mode from 
the full order system. In this way the error between 
the true system and the nominal system is nothing 
more than the second mode of the full order system. 
In addition, the real and imaginary parts of the first 
mode eigenvalue are both decreased by 0.10 to 
represent parameter errors in the nominal model. 

Model with Mixed Uncertainty 

This example addresses the ability of the model 
validation approach to generate an uncertainty model 
associated with the combination of parameter errors and 
unmodeled dynamics in the mass-spring-damper 
system. Figure 4 depicts the uncertainty model 
structure in block diagram form. Two uncertainty 
blocks are considered - parameter uncertainty and 
additive uncertainty. 

The additive uncertainty structure consists of (1) a 
full complex block A a in which each element of the 

2x2 matrix can take on any complex value such that 
the norm of the matrix is less than unity, and (2) a 
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single real scalar weight (i.e., = CO a l 2 ) that is 

chosen during the model validation solution. 

The parameter uncertainty is modeled as repeated real 
scalar perturbations. The weighting and uncertainty 
matrices have the following structure. 


“1 




Ap=l 


[1 °1 
to ij 


1 

0 


(02 


§ 2 


1 0 

Lo i Jj 

0 

1 °1 
0 


(43a) 


(43b) 


The first parameter 8 -| is a perturbation of the real 

part of the nominal model eigenvalue and the second 
parameter 82 is a perturbation of the imaginary part. 
Each Sj is real valued and can have arbitrary values 
such that |Ap|< 1. Each (Oj is also real but free to 


be chosen during the optimization. 

The solution of the model validation problem is 
accomplished using the optimization described 
previously over a range of frequencies. The input 
vector is [1 -0.5]' at each frequency. $ The objective 


function is the sum squared value of the uncertainty 
weights (Oj, i = 1 ,2,...,x . That is, the objective is to 


min 


|(01 

(0 T 



Figure 4 - Uncertainty model structure. 


: Laplace transform of 


1 

-0.5 j 


6(t) . 


For this example, the model error can be completely 
described by the chosen uncertainty structure. External 
disturbances and measurement noise are not considered 
and no weight is placed on the external disturbances. 

Figure 5 depicts results from the model validation 
solutions. The three optimal uncertainty weights (i.e., 
the additive uncertainty weight and the two real 
parameter uncertainty weights) are plotted at each 
frequency at which a solution was obtained. 

Figure 6 depicts the result of using the optimal 
uncertainty weights to generate a specific minimum 
norm validating uncertainty at each frequency point. 
The plot also shows several frequency responses of 
particular interest. The solid line is the frequency 
response of the "true system," the dashed line is the 
frequency response of the nominal model, and the 
dotted line is the frequency response associated with the 
error dynamics (i.e., the difference between the true 
system and the nominal model). 

The circles are the result of solving for minimum 
norm validating uncertainties and using them to 
compute the frequency response of the uncertainty 
model that corresponds to the error dynamics of the 
system, 

a{p 2 i[l-AP,i] 1 A Pj 2 1 • 



Figure 5 - Model validation uncertainty weights. 



Frequency (rad/sec) 

Figure 6 - Model validation solution. 
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The minimum norm validating uncertainties provide 
a good approximation of the error dynamics. Therefore, 
the model validation based uncertainty weights 
accurately characterize the differences between the "true 
system" and the nominal model. However, one might 
also expect the real parameter uncertainty weights to be 
0.1 at each frequency to reflect the actual parameter 
error. But this is not the case. So why are the real 
uncertainty weights so different from what might be 
expected? 

Note that the contribution of the parameter error to 
the error dynamics is most significant in the 
frequencies near the first mode while the contribution 
of the unmodeled dynamics is most significant at 
frequencies near the second mode. At frequencies near 
the first mode (where the parameter uncertainty 
dominates) the weights associated with the real 
parameter uncertainty are close to the theoretically 
exact value of 0.1. At frequencies away from the first 
mode the weights associated with the real parameter 
uncertainty are very small. 

The additive uncertainty is better able to characterize 
the uncertainty at frequencies away from the first mode 
than the parameter uncertainty. That is, the additive 
uncertainty permits validating signals for the system 
with a smaller weight than does the real parameter 
uncertainty. Recall that the solution approach seeks to 
find the smallest normed weights that achieve 
validation. The solution is not unique and a smaller 
norm solution than the theoretically exact solution was 
found to achieve model validation. 

Also note that because the parameter error is 
independent of frequency it should be the same constant 
value at all frequencies. Therefore, one way to improve 
the uncertainty model would be to solve the model 
validation problem again and force the parameter 
uncertainty weights to take on the maximum value 
obtained during the mixed uncertainty solutions. The 
sequential validation approach is an effective way to 
solve this problem. 

Sequential Model Validation 

The following example parallels the previous 
example. However, in this case it is assumed that the 
parameter uncertainty is known a priori. That is, the 
parameter uncertainty matrix corresponds exactly to the 
difference between the parameters in the true system 
and the nominal model. Recall that for the mass- 
spring-damper system there were two uncertain 
parameters, the real and imaginary parts of the first 
mode of the system, both of which were decreased in 
the nominal model by 0. 1 relative to their true values. 

The known parameter uncertainty was specified to 
be - 


0 . 1 | 


flpAp -| 


°-{o 


3 


(44) 


The additive uncertainty weights were then computed 
from the optimization of [3 to minimize the norm of 
the additive uncertainty weight and satisfy the 
constraint associated with its structure. 

Figure 7 depicts results from the sequential 
validation solutions. The three uncertainty weights 
(i.e., the optimal additive uncertainty weight and the 
two specified real parameter uncertainty weights) are 
plotted as a function of frequency. Note that, in the 
frequency range near the first mode, the additive 
uncertain weight is smaller than the corresponding 
uncertainty weight from the mixed uncertainty result 
(see Figure 5). This is expected because the larger 
parameter uncertainty is now better able to characterize 
the errors. 

The results of using the optimal uncertainty weights 
to compute a minimum norm validating uncertainties 
are shown in Figure 8. The plot shows the same 
frequency responses shown in Figure 6. However, the 
circles represent the uncertainty model dynamics 
associated with the optimized additive uncertainty 
weights with the real parameter weights specified. The 



Figure 7 - Sequential validation uncertainty weights. 



Frequency (rad/sec) 

Figure 8 - Sequential validation frequency responses. 
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minimum norm validating uncertainty provides a 
better approximation of the error dynamics than in the 
mixed uncertainty example (see Figure 6) as one 
might expect. 

Had the specified value of the parameter uncertainty 
been in error, the solution of the sequential validation 
problem could still have been accomplished. However, 
the resulting additive uncertainty would have had to 
describe the additional uncertainty in the assumed 
model due to the remaining parameter errors. In fact, 
this can be used as a means for selecting the best 
parameter uncertainty if it is not known a priori. 

The specified level of parameter uncertainty can be 
varied over a range of potential values. The sequential 
model validation method can then be used to compute 
the additional additive uncertainty required to achieve 
validation. A metric associated with the amount of 
additional uncertainty required to achieve validation can 
then be computed for each assumed level of parameter 
uncertainty. The parameter uncertainty producing the 
smallest additive uncertainty is the best choice for the 
mixed parameter / additive uncertainty structure. 

The results from an example of this approach for the 
mixed uncertainty case is shown in Figure 9. The 
metric for additive uncertainty is the area under the plot 
of the additive uncertainty weight for a range of 
specified parameter weights. The smallest area occurs 
for parameter uncertainty of 0.1 which corresponds to 
the least amount of additional additive uncertainty 
required for validation and is equal to the actual 
parameter error introduced into the nominal model. 

Concluding Remarks 

The solution of the model validation problem 
developed by Lim and Giesy provides an excellent 
foundation for generating uncertainty models for robust 
control design applications. The extensions to the 
Lim-Giesy approach developed herein provide an 
ability to address some additional aspects of the 
uncertainty modeling problem. The extended 
optimization framework allows the determination of 
the best (in a normed sense) combination of uncertain 
elements (e.g., parameter and additive uncertainty) to 
achieve validation. 

The sequential solution approach provides a means 
by which uncertainty structure can be addressed in the 
context of the standard model validation solution 
framework. This approach provides more control over 
the way in which various uncertain elements can be 
combined and does not require the engineer to rely 
solely on the optimizer to find the 'best' solution. 

The sequential solution approach has some added 
benefit when there is more knowledge about the true 
system than just input/output data. For example, the 



Figure 9 -Parameter uncertainty selection with 
sequential validation. 

situation when certain parameter variations are 
known. This knowledge can be included in the model 
validation solution so that the additional uncertainty 
model components will only have to characterize 
those aspects of the system not otherwise described 
by the known components. The result is less 
conservative uncertainty models and thereby better 
controller designs. 
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